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Abstract. Our report are the 2.5D and 3D numerical models of the interaction 
of the solar wind with the Solar Probe Plus (SPP) spacecraft. These results should 
be interpreted as an engineering document for which the derived SW interaction with 
SP+ could have consequences for both plasma wave and electron plasma observations. 
The SPP model includes 3 main parts, namely, a non-conducting heat shield, a support 
system, and cylindrical section or spacecraft bus that contains the particle analysis 
devices and antenna. The simulation was performed using 3D hybrid code which 
describes proton dynamics by particle motions whereas the electrons are considered in 
a fluid approximation. The results of the quasi-parallel and quasi-perpendicular solar 
wind - SPP interaction at the distance r = 4.5 R s were discussed in the preliminary 
Report (Lipatov, Hartle and Sittler, 2009). In present Report we discuss the results 
of simulation of the solar wind - spacecraft quasi-parallel interaction at the distance 
r = 9.5 R s . In this regime the value of the ion thermal velocity is about of the value 
of the bulk velocity and the polarization electric hied is strong enough so that only 
a very small cavern is formed behind the heat shield and the bus. One observes an 
excitation of the low frequency Alfven and whistler type wave directed by the magnetic 
held with an amplitude of the electromagnetic held oscillation about of (0.06-1.5) V/m. 
The compression waves and the jumps in an electric held with an amplitude of about 
(0.15-1.5) V/m were also observed. The results of our hybrid simulation will be useful 
for understanding the plasma environment near the SP+ spacecraft. Future simulation 
will take into account charging of the spacecraft, charge separation effects, outgassing 
from the heat shield, photoionization and electron impact ionization effects near 
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the spacecraft. We also need to perform a longer simulation in order to receive a 
’’steady-state” solution with formation of a plasma wake in a quasi-parallel interaction. 
We are now in a position to do this, since these codes have now been converted to 
the hybrid MPI/OpenMP environment and can now run on the Pleiades and Discover 
massively-parallel supercomputers with distributed memory. 

Numerical simulation studies; Key words: Solar wind; Spacecraft; Alfven waves, 
Whistlers, Ionospheres; Atmospheres; Induced magnetospheres; Magnetic barrier; 
Satellites 
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1. Introduction 

The numerical simulation of the interaction of space plasma with spacecraft devices 
plays a key-role in their design. The study of the excitation of waves near a rapidly 
moving body has long history starting from the quasi-analytical approach (see e.g. 
Alpert (1974) and the references thereby) to the recent electrostatic modeling. While 
the electrostatic studies produced very important information concerning the charging 
of spacecraft, plasma void and wake, the problem of excitation of the electromagnetic 
waves upstream and downstream flow near spacecraft are great interest for onboard 
measurement in the SW, magnetosphere and ionosphere. This is especially true for 
SPP where new regimes of solar wind plasma are being explored. The simulation 
provides the general characteristics of the plasma environment and the electromagnetic 
field distribution near the probe and at it’s surface. The external surface of the Solar 
Probe Plus consists of a conducting part and an insulating part. It moves in the 
supersonic/subsonic and superalfvenic/subalfvenic solar wind flow. 

In the first approximation the plasma environment near the Solar Probe+ may be 
similar to the plasma environment near the Moon with a weakly conducting surface. In 
the case of a non-conducting model of the Moon the solar wind particles penetrate the 
surface on the day-side of the Moon whereas on the night-side a plasma wake with the 
low-density void is formed. On the day side of the Moon’s surface the plasma particles 
are absorbed, and the perturbation region forms a thin boundary layer of thickness 
S ~ c/a;p e oc 1km at 1AU. (Neugebauer, 1960). On the night side the perturbation 
region is bounded by a surface of weak perturbations forming a cone with half-apex 
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angle sin 2 9 = 1 / M\ + 1/M| (Ma oc 8 and Ms oc 8 are the Alfven and sound Mach 
numbers, Va ~ 50km/s, w t h ~ 50km/s, [/ ~ 400km/s, and d ~ 10°) (Whang, 1969; 
Lipatov, 1976). 

No bow shock or oblique Mach cone are formed. The perturbation of the magnetic 
held inside the plasma wake is mostly determined by diamagnetic current and the drift 
currents due to a curvature of the magnetic held and acceleration of the bulk velocity 
[Wang, 1968; Wang, 1969; Wang and Ness, 1970; Lipatov 1976, Lipatov, 2002], Kinetic 
instabilities may also play a role in wake dynamics. Hybrid simulation of the interaction 
between supersonic plasma how and a weakly conducting body shows formation of a 
strong oblique double shock structure inside the plasma wake [Lipatov, Motschmann, et 
ah, 2005], 

2. Formulation of the Problem and Mathematical Model 

2.1. Solar Wind Simulation Model 

The interaction of solar wind particles with the ” Solar Probe+” is more complicated 
because the gyroradius of the protons (for Maxwellian core of the velocity distribution 
50 — 500 km) is much larger than the size of the spacecraft (2-3m). For energetic 
components the ion gyroradius may be extremely large with respect to the spacecraft. 
So the kinetic approximation for ions is an essential part of the mathematical model. 
The electron gyroradius may vary in range from 2.5 x 10 2 m at Earth orbit to much 
smaller scale, 1 m, near Sun. Therefore, electrons need a huid-kinetic approximation to 
study the plasma environment along the ’’Solar Probe+” trajectory. 
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We study the plasma environment near the ” Solar Probe+” by various simulations; 
e.g., (a) standard hybrid simulation (ion in kinetic approximation, electron in fluid 
approach) on the large scale; (b) fully kinetic implicit simulation with kinetic model 
for electrons and ions (see e.g. Hewett and Langdon, 1987; Lipatov and Lobachev, 
1996; Lipatov, 2002; Damiano, Sydora, Samson, 2003; Lipatov and Rankin, 2008) 
incorporated in the large scale hybrid model. The last simulation will take into account 
charge separation near the surface of the spacecraft and finite electron gyroradins 
effects. We will take into account the realistic distribution of the spacecraft surface’s 
conductivity. Our simulation will serve as an expert system for design of the ’’Solar 
Probe+” spacecraft. The present model of the interaction of the solar wind with 
the SP+ does not take into account several effects in plasma environment near the 
spacecraft. Future simulation will take into account the charging of the spacecraft, 
charge separation effects, ontgassing from heat shield, photoionization and electron 
impact ionization effects near the spacecraft. 

To study the interaction of the solar wind with the SP+ spacecraft we use a 
quasi-neutral hybrid model for ions and electrons. 

In our coordinate system the X (radial) axis is parallel to the solar wind velocity - 
U 0 , Z is aligned with the equatorial plane, and Y = Zx X, Fig. 1. 

In the hybrid simulations described here, the dynamics of upstream ions is 
described by kinetic approach, while the dynamics of the electrons is described by a 
hydrodynamical approximation. 

The single particle ion distribution function /j(t, x, v) has to satisfy the 
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Vlasov/Boltzmann equation 


— f-L —f _l F 19 

0T* + V dx + 



( 1 ) 


where F symbolizes forces due to electric and magnetic fields acting on the ions, F co u 
is the collision term. In this paper we use the particle-mesh model for ion dynamics 
instead of the Vlasov equation, Eq. ( 1). 

The single ion particle motion is described by the equations (see, e.g. Eqs. (1) and 
(14) from [Mankofsky, Sudan and Denavit, 1987]): 


dr it i dvjj e / v,,i x B \ m e u ie 

dt 1,1 ) dt c ) Mieni 


( 2 ) 


Here we assume that the charge state is Z t — 1 and that all ions have the same mass 
Mi. U i and J denote the charge-averaged velocity of all (incoming and pickup) ions and 
the total current, Eq. (6). The subscript i denotes the ion population and the index l 
is the macro-particle index. u ie is collision frequency between ions and electrons, that 
may include Coulomb collisions and collisions due to particle-wave interaction. Note 
that the collision rates used in Eq. (2) must depend on the individual velocities of ions 
and electrons. However, we use the effective resistivity rj, r) = a ^ 1 = m e /(ne 2 T e ), where 
T e = z/" 1 . The electrical conductivities may be estimated as 


a± = a{Tl /2 , ct|| = 1.92oj_, oi = 0.9 x 10 13 /((A/10)Z;) s" 1 • eV" 3/2 , (3) 


where T e denotes the electron temperature and A is the Coulomb logarithm (see, e.g., 
pages 215-216 from [Braginskii, 1965]). For typical solar wind parameters T e = 100 eV 
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(electron temperature), no = 10 4 cm -3 (density), and electrical conductivities are 
a _ l ~ 4.7 x 10 13 s _1 and a\\ ~ 9.2 x 10 13 s _1 . 

In our simulation we use the value of the effective conductivity that is much smaller 
than the real value to suppress ’’shot” noise. The numerical ’’shot” noise is a result of 
a fluctuation in density and bulk velocity due to a small number of particles per cell in 
the particle in cell simulation. We also use an effective conductivity for modeling Solar 
Probe+’s bus. Hence, we must keep the first collision term in the right hand side of Eq. 
( 2 ). 

In the nonradiative limit Ampere’s law is given by 

47 r 

— J = V x B; (4) 

c 

and the induction equation (Faraday’s law) by 

+ VxE = 0. (5) 

c at 

The total current is given by 

J — J e T Jj , J; — ( !) j U, . (6) 

We further assume quasi- neutrality 

n e « Hi. (7) 


The equation of motion of the electron fluid takes the form of the standard 
generalized Ohm’s law (e.g. Braginskii, 1965): 


1 / , i 1 „ m e J 

E = (J c x B) Vp e H W* — 

en e c en e e ne 


m d 
e dt 


U e , 


( 8 ) 
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where p e = nm e (v ®) / 3 = n e /csT e , and v' e are the scalar electron pressure and the thermal 
velocity of electrons, and the electron current is estimated from Eq. (6). The last term 
on the right side of Eq. (8) is the electron inertial term. At 5 R s we have no = lCrcm -3 , 
v’ e = 5000 km/s, and from Eq. (3) one gets cr e>i = 0.4 x 10 14 s -1 for background plasma 
parameters. 

Since we suppose that electron heating due to collisions with ions is very small, the 
electrons are considered as adiabatic fluid. For simplicity we assume that the electron 
pressure may be represented as: 

p e oc nf /3 . (9) 

We also assume that n e ~ n\. Otherwise, we have to calculate the electron pressure from 
heat balance for electrons (see, e.g., Braginski, 1965) taking into account the heat fluxes. 
The hybrid model allows us to study the dynamics of the ion velocity distribution and 
the effects of a finite ion gyroradius. The electron gyroradius and charge separation 
effects are not resolved in hybrid models and we have to use a full kinetic model in 
future simulations. 

2.2. Solar Probe Plus Spacecraft Model 

Figure 1 shows a scheme of the Solar Probe Plus (SP+) spacecraft, a system of 
coordinates and a direction of the plasma flow and magnetic field. The SP+ model 
includes 3 main parts, namely, a non-conducting heat shield (1), a support system (2), 
and a cylindrical section or spacecraft bus (3) that contains the particle analysis devices 
and antenna. The heat shield has the following geometrical parameters: diameter, 
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D shield = 2.7 m, and thickness = 0.335 m. The gap between the heat shield and 
cylindrical section is 1.188 m. Note the gap is now covered with radiators, so there is 
no penetration of the solar wind through the radiators. The cylindrical section or bus 
has a diameter of 1.026 m and length 1.188 m. We also take into account the effective 
resistivity of the Solar Probe+’s bus, the heat shield and the gap between the heat 
shield and the spacecraft bus: 

Pbus = (1 — 15) x 1(T 3 ohm-m ; Pshieid = (3 — 15) ohm-m ; p gap = (3 — 15) x 1CT 3 ohm-m 

( 10 ) 

Our code solves equations (1) - (5), (6) - (9) (Lipatov, Hartle, Sittlcr, 2009). Here we 
note that the gap has several trusses to provide a mechanical interface between bus and 
heat shield. If electrical conductivity of trusses is high enough that differential charging 
between spacecraft bus and heat shield is low, then the heat conduction from heat shield 
to spacecraft bus may be too high. We still don’t know the exact value at this time with 
regard to truss electrical conductivity parameters. 

Initially the computational domain contains only supersonic or subsonic solar 
wind flow with a homogeneous spatial distribution. The magnetic field and electric 
holds are B = Bo and E = Eo = — Uo x Bo. Inside the Solar Probe+ spacecraft the 
electromagnetic holds are E =0 and B = B 0 , and bulk velocities of ions and electrons 
are also equal to zero. In the cases examined here we choose for the magnetic held a 
spiral angle 0b u = 11°. 

Far upstream (x = —DXj 2), the ion hux is assumed to have a Maxwellian 
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distribution, 

/ = “ 3/2 exp ~ ^ V 2 U ^ , (11) 

zv th 

where Vth and U are the thermal and the bulk velocities of the solar wind plasma flow. 
We have not included a contribution from the field aligned strahl electrons which are 
non-Maxwellian. 

Far downstream, we use a buffer zone to provide the return of particles with 
negative values of the u velocity component. We also adopted Sommerfeld’s radiation 
condition for the magnetic held. On the side boundaries (y = ±DY/2 and z = ±DZ / 2), 
unperturbed boundary conditions were imposed for the incoming flow particles and the 
electromagnetic field. At the spacecraft surface the particles are absorbed. There is no 
boundary condition for electromagnetic fields, and we also use our equations for the 
electromagnetic field inside the spacecraft but with internal conductivity and the bulk 
velocity that is calculated from the particles. In this way the jump in the electric held 
is due to the variation of the values of the conductivity and bulk velocity across the 
surface of the spacecraft PS+. The position of the center of the bottom of the heat 
shield is defined asa: = 0,|/ = 0,2; = 0. 

The three-dimensional computational domain has dimensions DX = 10 L, DY = 8 L, 
and DZ = 8 L, where L = D shie i d / 2 = 1.35 m is the radius of the heat shield. We use a 
mesh of 401 x 301 x 301 grid points, and 9 x 10 8 particles for protons in a homogeneous 
mesh computation. We use a subcycling procedure for time integration of the particle 
and the electromagnetic equations (see, e.g. Lipatov 2002). The time step for the particle 
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update, At p satisfies the condition v max A t p < mm(Ax, Ay, Az)/16, whereas the time 
step for the electromagnetic held time integration, Atf, is much smaller than the time 
step for particle update. Atf satisfies the condition v max Atf < min(Ax, Ay, Az)/6400. 
The grid spacings are as follows: Ax = Ay = Az = 0.036 m. 

The relationships between dimensional (U, E, B, p e , n, T) and dimensionless (TJ 7 , 
E', B', p' e , n' , T') parameters may be expressed via the dimensional upstream values as 
follows: 

U = XJ'Uq, E = E'BqUo/c, B = B'B 0 , p e = Pe'PeO, 

n = n'no, T = T'M-JJq, (12) 

whereas the dimensional time and distance may be expressed via the bulk velocity Uq 
and characteristic scale L = D s h ie id/ 2: 

t = t'L/Uo, x = x'L. (13) 

The global physics in SP+’s environment is controlled by a set of dimensionless 
independent parameters such as Alfven Mach number M\, ion and electron plasma betas 
Pi, P e , electron/proton mass ratio rri/ M p , diffusion lengths, and the gyroradins e = p C i/L. 
Here p ci = U 0 /(eB /Mjc) = M\c/cu pi and the ion plasma frequency ca pi = yj Airnoe 2 / M-,. 
The actual value of the proton gyroradins is about (0.31 — 2.5) km using the above 
formulas. The grid spacing has the value A x = L / 50 = 0.036 m. In order to model the 
ion kinetics, the ion gyroradins must be resolved on the grid. 
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2.3. Numerical Method 

We employed a standard particle-in-cell (PIC) method with a homogeneous grid. 
The time integration of the particle equations of motion uses a trapezoidal leapfrog 
scheme (see e.g., Lipatov, 2002). The time integration of the electromagnetic field 
equations uses an implicit finite difference scheme (see, e.g., Lipatov (2002)). 

Equation (8) is evaluated at a time level between the n and n + 1 levels. The 
electric and magnetic fields at this time level, n + 6 , are given by 


E a+e = eE n+l + ^ _ 0) E n 


(14) 


B n+e = 9B n+1 + (1 - 0)B n . 


(15) 


In the predictor step one needs to calculate E and B at time level n + 6, and for 
this purpose we used second Maxwell equation, which gives 


H“'" = B” - ScAiV x E n+0 . 


(16) 


The finite difference approximation for the electric field at time level n + 9 in 
combination with (16) and the electron velocity equations may be produced by the 
following way: 

AE n+e + (V x B n+e ) x I 


+S((V x (V x E n+0 )) x B n ) = Q, 


(17) 


where 


9 = 


9 Ate 

Mp 


d = 9n e l^At, A = n e , 


(18) 
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and 


N s 


Q = 


ej: 

vfe=l 


e/3, 


k = 1 


Mi 


V x B n x B' 


2 Mi 


Vp e + n e l^V x B' 


U, = Ui - eJ/(M>,). 


(19) 

( 20 ) 


( 21 ) 


Note that we dropped the electron inertia term in the above finite-difference scheme for 
the electric field. 

In the case with adiabatic electrons one can split the total electric field into the 
sum of inductive ( E \ ) and electrostatic (E 2 ) fields: 


E — Ei + E 2 , 


( 22 ) 


where E 2 satisfies the condition 

Vx£ 2 = 0. (23) 

Then we can solve (17) for component E\ neglecting Vp e - The electrostatic (polarization) 
electric field E 2 can be calculated from p e (Lipatov, 2002) because 

E 2 = — — Vp e . (24) 

en e 

Let us consider two meshes. The first mesh contains the nodes which are 

located in the center of a cell. It is used for the computation of the density, current, 
bulk velocity, electron pressure and inductive electric field. The second mesh contains 
the nodes i ± 1/2, j ± 1/2, k ± 1/2, which are located at the corners of a cell. This 
mesh is used for the computation of the magnetic field components, electrostatic field 
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and final electric field. Let us assume that Jf +1,/2 , o; pe and n ; are known at time level 
n+ 1/2. Then, expressing V x E'i +e and V x (V x E™ + ) via central finite differences 
at each cell center, one can obtain the following 3x3 matrix equation: 

L ■ = F ' <*) 

Equation (25) may be solved by iteration. In each iteration the electric field on the right 
side is given. The iteration continues until some convergence criterion is satisfied. At 
the same time the electrostatic field E 2 is calculated at from Eq. 24. 

In the corrector step, when E 2 is obtain, one can update the magnetic field, using 
again the second Maxwell equation: 

nAf 

B n+1 = gn _ __y x ( E n+1 + (26) 

The second term on right side of Eq. 20 keeps a large parameter, so that a strong 
shot noise in a density computation may cause a strong oscillation in the electric field. 

Since the gyroradius must be resolved, a grid point spacing of less than 1 gyroradius 
is required in order to avoid numerical dispersion and dissipation. On the other hand, 
good statistics are required, therefore a sufficiently large number of particles per cell 
have to be used (i.e., to obtain low “shot” noise, which manifests itself as fluctuations 
in the numerical plasma parameters due to a small number of particles per cell). 

A multiscale simulation will be conducted using adaptive mesh and particle 
refinement techniques with composite grids. This code exploits a splitting-of-particles 
procedure to keep a low level of “shot” noise on the hirer mesh. At the end of the global 
time step we use the synchronization for the values of the electromagnetic held on coarse 
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and finer meshes. A.S. Lipatov has developed a version of the Complex Particle Kinetic 
method (see, Hewett, 2003) [in collaboration with D.W. Hewett (LLNL) and M.R. 
Combi (Univ. of Michigan)] which allows us a reduction in computational resources 
by factors of 100 - 1000 in comparison with standard PIC simulation (Lipatov, 2008a). 
These codes were optimized for parallel computation using MPI and OMP. 

One 3D hybrid simulation takes around five months for computing with 24 
processors and 64-128 GB RAM on the ” shared memory” Columbia SGI system 
in the NASA Ames Supercomputer Center. Currently, our code operates in the 
parallel processing with the OpenMP environment. Converting this code to the MPI 
environment was performed by technical staff members of the Information Technology 
and Visualization Office (NASA GSFC) starting from Spring 2008 and we are now testing 
the MPI and OpenMP versions of the code on Pleiades and Discover massively-parallcl 
supercomputers with distributed memory. 

3. Results of the Simulation 

To study the interaction of solar wind with the SP+ we use the following sets of 
solar wind plasma parameters: — 0.1 (T t = IMA'); /3 e = 0.1 ( T e = IMA'); magnetic 

held, B 0 = 1500 nT; bulk velocity U 0 = 200km/s, density nsw = 5 x 10 3 cm -3 , Alfven 
Mach number Ma = 1. — 1.5 and 9b u = 11°. The value of the plasma bulk velocity in 
the spacecraft frame may vary from 200km/s to 800km/s so we can use a higher Alfven 
and sonic Mach numbers in our simulation. Table 2 shows the all parameters that 
correspond the cases: (a), (b), (c), (d), (d’), (e) and (e’). These parameters correspond 
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to the values of the motional electric field Eq ps 0.3 V/m for r = 9.5 R s . For the distance 
from the Sun of about 5 R s the value of the magnetic field is of about 6000 nT and the 
solar wind density is of about 2 x 10 4 cm -3 (Lipatov, Hartlc and Sittler, 2009). We 
present here the results of simulation for the distance, r = 9.5 R s . The time step for 
electromagnetic field update in the simulation of the quasi-parallel interaction is 48 
times smaller than the time step in an oblique case. In the present simulation, the 
dimensionless scale of proton gyroradius is e = p p /L = 10 3 . We discuss the results of 
3D hybrid simulations in Sect. 3.1. The results of 2.5D simulation will be discussed in 
Sect. 3.2. 2.5D simulations include both the particle absorption boundary condition 
at the heat shield (case d) and the particle reflection boundary condition on the heat 
shield (case e). We are checking the results of 2.5D simulations in cases (d) and (e) with 
simulation having smaller grid spacing (case d’ and case e’) but with the same plasma 
parameters and boundary conditions. 

3.1. Quasi-parallel Interaction of the Solar Wind with Solar Probe Plus. 3D 
Simulations 

In this section we discuss the results of 3D simulations for three cases, namely, 
case (a): lower number of the macro-particles, N p « 2.7 x 10 9 , so there are about 220 
macro-particles per cell; (case b): a large number of the macro-particles, N p & 8 x 10 9 , 
so there are about 700 macro-particles per cell with extra-smoothing of the polarization 
electric field; (case c): same as case (b) except now without extra-smoothing of the 
polarization electric field. It seems that the larger number of the macro-particles is 
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not large enough to produce a solution with a low ’’shot” noise for plasma parameters 
near r = 9.5 Rs- Note, that a ’shot” noise is a statistical effect of a fluctuation of the 
plasma parameters such as bulk velocity, density etc. which results in the fluctuation 
in the electromagnetic field. These fluctuations in the electromagnetic held produce 
overheating of the macro-particles. For that reason we use a special smoothing procedure 
to reduce this noise (case b). 

3.1.1. Linear Regime. The results for a long time simulation show the 
formation of only a small cavity with low plasma density. The reason is the large 
thermal speed in comparison with the bulk velocity and a strong polarization electric 
held, E po i = —Vp e /en e , in cases (a), (b), (c) with j3 e = 0.1. A large polarization electric 
held can support a strong gradient in the density prohle and the result is the formation 
of a large electric held behind the bus and the heat shield. 

Figures 2 and 3 show 2D cuts for proton density in the y — x and z — x planes 
at t = 0.014T ce ( t = 0.0067t rans ;t), where T ce is the electron gyroperiod and T tran sit 
denotes the time for particle transit from the left boundary to the right boundary of the 
computational domain, (cases a and b). We see the picture of the solar wind how around 
the spacecraft. Behind the heat shield a density prohle forms a cone due to expansion 
of the external plasma into the gap between the heat shield and the spacecraft bus. 
(Present SP+ spacecraft from AO has this region closed oh with radiators). Behind the 
spacecraft bus there is no apparent plasma wake due to low value of the Alfven Mach 
number, Ma = 1.5, and relatively high plasma beta values, f3 p = 0.1, /3 e = 0.1 (Table 2). 

In these simulations the current closure near the spacecraft is very complicated 
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and is directed by the external magnetic field. Figures 4-9 show the 2D cuts for the 
electric field at t — 0.014 T ce . A linear perturbation occurs in the electric field in the 
form of whistler/Alfven waves in upstream and downstream regions at the beginning of 
the simulation. 

In case a, the value of the perturbation in the electric field is about 6E & 

(0.02 — 0.6) Eq, for which Eq « 0.3 V/m is defined as a standard unit for these 
computations (see Table 3 for a summary of the electric and magnetic field perturbation 
amplitudes) . 

Figures 10 - 15 show 1-D cuts for the electric field profile along the x axis through 
the point at y — 0 and z = 0, y axis through the point at x — 1.5 L and z — 0, and z 
axis through the point at x — 1.5 L and y — 0. 

In case a, shown in Fig. 10, far from the spacecraft, the value of perturbation in 
the electric field component, E x , is about of 8E X ~ 0.01 E 0 , whereas 5E y m 0.05 E 0 and 
SE Z ~ 0.05 Eq, in the x- direction (upstream and downstream). We assume that these 
waves represent the right-polarized whistler with wavelength A ~ 2 m. We also observe 
this type of waves in the downstream region behind the bus. 

In y- direction, Fig. 12, the perturbation in the electric field far from spacecraft are 
the following: E x « 0.003 E 0 , E y & 0.01 E 0 and E z « 0.005 Eq. 

We also observed the excitation of a high frequency wave that propagates across 
the magnetic field in the z - direction, Fig. 14, with wavelength A ~ 0.15 m. The 
perturbations in the electric field are the following: E x « 0.04 Eq, E y « 0.05 E 0 and 


E z m 0.2 Eq 
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In case a, near the spacecraft the value of perturbation in the electric field 
components, E y and E~, is about (0.05, 0.15) E 0 (whistler type waves) in the x- direction 
(upstream and downstream) in Fig. 10; 5E X ~ 0.02 E 0 , 5E y ~ 0.4 E 0 and 5E Z ~ 0.1 E 0 
in the y- direction (Fig. 12); and SE X ~ 0.06 E 0 , 5E y ~ 1.0 E 0 and 5E Z ~ 0.4 E 0 in the 
2 - direction (Fig. 14). Near the surface of the spacecraft, the jumps in the electric field, 
5E X = (—0.075, —0.15) E 0 , 5E y = (-0.075, 0.02) £ 0 and 5E Z = (-0.12, 0.15) £ 0 occur 
behind the heat shield (Fig. 10). Behind the bus (Fig. 12), the jumps are SE X = 0.01 E 0 , 
5E y = —0.2 E 0l 5E Z = — 0.025 E 0 . At the side surface of the bus (Fig. 14), 

5E X = (—0.005, 0.025) E 0 , 5E y = (—0.1, —0.14) E 0 and 5E Z = (—0.07, 0.01) E 0 in the y - 
direction and 5E X = (—0.02, 0.03) E 0l 5E y = (—0.3, —0.4) E 0 and 5E Z = (0.1, —0.175) E 0 
were observed in the 2 - direction. 

In Figures 11, 13 and 15 we show electric field for case (b). 

Far from the spacecraft, the value of perturbation in the electric field component, 
E x , is about of 5E X « 0.1 E 0 , whereas 5E y & 0.2 E 0 and SE Z ~ 0.3 E 0 , in the x- direction 
(upstream). We assume that these waves represent the right-polarized whistler with 
wavelength A ~ 0.2 m. 

In y- direction, the perturbation in the electric field far from spacecraft are the 
following: E x ~ 0.005 E 0l E y ~ — 0.2 E 0 and E z ~ 0.01 E 0 . 

In 2 - direction, the perturbations in the distant electric field are the following: 

E x « 0.05 E 0l E y ~ —0.2 E 0 and E z « 0.02 E 0 . 

Near the spacecraft, Fig. 11, the E x profile has a jump, 5E X = — 0.1 -Fq before the 


heat shield and 5E X = —0.3Eq before the bus and 5E X = 0.1 Eq behind the bus. Near 
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the spacecraft we see an excitation of a weak oblique Alfven (magnetosonic) waves with 
an amplitude, 8E X « ±0.05 E 0 along y axis (Fig. 13), and 8E X r* — 0.05-F 0 along the z 
axis (Fig. 15). The prohle of the E y component has a jumps, 5E y ~ ±1 Eq in y direction 
(Fig. 13) and SE y r; 0.4 E 0 in z direction, (Fig. 15). The components E z has a jump, 
8E Z ~ ±0.3 Eq in the y (Fig. 13) and 5E Z ~ ±1 E 0 in the z directions (Fig. 15). 

Figures 16 - 21 demonstrate 2D cuts for the magnetic held. At the front of the 
heat shield a formation of the magnetic held barrier is observed. At the side parts 
of the computational domain the whistler/ Alfven waves were observed in the x — z 
plane. The value of the perturbation in the magnetic held component, B y . is about of 
5B y ~ 0.005 B 0 (« 7.5 nT for B 0 « 1500 nT.). The value of the perturbation in the 
magnetic held B z is about of SB rs 0.02 Bq. 

3.1.2. Regime of Nonlinear Saturation. Let us consider the results of 
simulations at later times, t = 0.1 T ce (case a), t = 0.15T ce (case b), t = 0.29 T ce (case c). 

Figures 22 - 24 show a 2D cut for proton density in the y — x and z — x planes. The 
small voids with low plasma density are observed behind the heat shield and the bus. 
However, the size of these is limited by a strong polarization held. Possibly we need a 
several times longer simulation. 

Figures 25 - 33 show the 2D cuts for the electric held. Since the density distribution 
had not changed during the simulation and the perturbation in the electromagnetic 
held reaches the saturation level, we do not find any strong transverse perturbation 
in the electric helds in the upstream region, but we do see large perturbations in the 


downstream region. 
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Figures 34 - 45 show 1-D cuts for the electric held profile along the x axis at the 
point y — 0 and z = 0, along the y axis at the point x = 1.5 L and z — 0 and along the 
z axis at the point x = 1.5 L and y — 0. 

In case a, Figs. 34 and 37, the value of these perturbations far from the spacecraft 
is about of SE r; (3.0 — 5.0) E 0 r (7.0 — 10.0) V/m in x - direction (upstream). 

Far from the spacecraft in the downstream region, the value of perturbation in the 
electric held component, 8E X , is about of SE X ~ ±(3.0 — 5.0) E 0 in the x - direction, 
SE X rs ±(3.0 — 5.0) E 0 in the y - direction (Fig. 40), whereas SE X ~ ±(3.0 — 4.0) E 0 in 
the z- direction (Fig. 43). The value of perturbation in the electric held components, 
E y and E z , is about of 5 E r ±(3.0 — 5.0) E 0 (whistler type waves) in the x- direction 
(downstream), Figs. 34 and 37. SE y ~ ±(3.0 — 4.0) E 0 in the y- direction (Fig. 40), and 
SE y r ±2.0 Eq in the z- direction, Fig. 43. The value of perturbation in the electric 
held component, E z , is about of 8E Z r ±(3.0 — 5.0) E 0 in the y- direction, Fig. 40, and 
8E Z r ±(3.0 — 5.0) E 0 in the z- direction, Fig. 43. Near the spacecraft, Fig. 43, the E x 
prohle has a jump, SE X = —250 E 0 behind the cylindrical section (conducting bus). We 
also see a jump in E y component of the electric held, 8E y R ±55.0 Eq in y direction 
(Fig. 40) and SE y r (7.0 — 8.0) E 0 in z direction. The component E z has a jump, 

8E Z r; (—60.0, —400.0) E 0 in x direction (Fig. 34), 8E Z & (4.0, — 6.0) E 0 in y direction 
(Fig. 40) and 8E Z ~ ±60.0 E 0 z direction (Fig. 43). 

In case b, far from the spacecraft, the value of perturbation in the electric held 
component, E x , is about of 8E X ~ (1.0, 0.5)£J 0 i n H ie x ~ direction (upstream and 
downstream), Figs. 35, 38, whereas 8E X ps ±0.5 Eq in the y- direction, Fig. 41, and 
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5E X « (1.0, 0.5) Eq in z - direction, Fig. 44. The value of variation in the electric 
field components, E y and E~, is about of ±0.5 E 0 in x, y and z (upstream), Fig. 38. 
Downstream, the value of the perturbations in E y and E z are about of SE « 2.0 Eq in x 
- direction, Fig. 38. 

Near spacecraft, Fig. 35, the E x profile has a jump, 5E X ~ — 190 A 0 in x - direction 
behind the heat shield and the bus, and 8E X & (0.5 — 1.5) E 0 along y (Fig. 41) and z - 
axes (Fig. 44). We also see a jump in E y component of the electric field, SE y ~ 100 E 0 
in x - direction behind the bus (Fig. 35), SE y « ±40 Ao in y - direction (Fig. 41) 
and SE y « (6.0, 7.0) E 0 in z - direction (Fig. 44). The component E z has jumps, 

SE ~ — 160 E 0 in x - direction (Fig. 35), SE Z ~ — (5.0, 6.0) E 0 in y (Fig. 41), and 
5E Z « ±40 A 0 and z - directions (Fig. 44). 

In case c, far from the spacecraft, the value of perturbation in the electric 
field component, E x , is about of SE X ps ±2.0 E 0 in the x - direction (upstream and 
downstream), Fig. 39. Similar oscillations occur for cuts along y - direction (Fig. 42) 
and z - direction (Fig. 45). The value of perturbation in the electric field components, 
E y , is about of 5E y « ±2.0 Eq ( whistler/ Alfven waves) in the x- direction (upstream 
and downstream), y, and z- directions, Figs. 39, 42 and 45, respectively. The value 
of perturbation in the electric field component, E~, is about of 5E Z ps ±2.0 E 0 in the 
x, y and z - directions, Figs. 39, 42 and 45. Near spacecraft, Fig. 36, the E x profile 
has a jump, SE X ps — 230 E 0 in x - direction behind the heat shield and the bus. 

Near spacecraft we see an excitation of a strong Alfven waves with an amplitude, 

SE X ~ (2.0, —3.0) Eq along y - axis (Fig. 42), and SE X ~ (—1.0, —2.0) Eq along the z - 
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axis (Fig. 45). We also see a jump in E y component of the electric field, 5E y « — 250 E 0 
in x - direction behind the bus and SE y ~ ±40.0 E 0 in y - direction and SE y ~ 4.0 E 0 
in z - direction, Figs. 39, 42 and 45, respectively. The component E z has jumps, 

SE « — 100.0 Eq in x - direction behind the bus (Fig. 36), SE « — 4.0-Eo in y - direction 
(Fig. 42), and SE « ±35.0 E 0 in z - direction (Fig. 45). As before, the above figures 
demonstrate a strong oscillation in the electric held. We deducted the polarization 
electric held from the total electric held and we found the regular current closure that 
corresponds the Alfven wing. 

The above figures demonstrate a strong oscillation in the electric held. We also 
deducted in case (a) the polarization electric held from the total electric held and we 
found the regular current closure that corresponds the Alfven wing. If we assume the 
’’shot” noise Sn/n = 0.07 then a fluctuation in the polarization electric held may be 
estimated as SE po i = —Wp e /en e k 2 x 10 2 E 0 . The nature of the ’’shot” noise is a 
statistical fluctuation of the particle velocity distribution in cell due to a limited number 
of the macro-particles per cell (see e.g., Birdsall, Langdon and Okuda, 1970). This value 
represents the upper limit for the polarization electric held, however, the level of electric 
held fluctuations is much lower because we use a ’’quiet start” in our simulation i.e. we 
started from a plasma distribution without any fluctuation in the density and current. 
In case (b) the ’’shot” noise in density is estimated as Sn/n = 0.04, then a fluctuation 
in the polarization electric held may be estimated as SE po i = — ' Vp e /en e ~ 10 2 E 0 . 
And, finally, in case (c) the estimation of the ’’shot” noise gives the same value in the 
polarization electric held as in case (b). 
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Figures 46 - 54 demonstrate 2D cuts for the magnetic field. In case (a), we do 
not find any magnetic field barrier near the heat shield. At the side parts of the 
computational domain, compression waves occur in the simulations. The value of 
perturbation in the magnetic field is about of SB « 0.02 Bq. In case (b), at the front 
of the heat shield a formation of the magnetic field barrier or build up was observed. 
The formation of the Alfven wing was also observed. The value of perturbation in the 
magnetic field is about of SB ~ (0.004, 0.03) B 0 . In case (c), at the front of the heat 
shield a formation of the magnetic field barrier or build up was observed. At the side 
parts of the computational domain, compression waves were observed in simulations. 
The value of perturbation in the magnetic field is about of SB « 0.02 Bo- 

Summary 1. 3D hybrid simulations show the following global plasma dynamics 
near the SP+ spacecraft. The heat shield and bus generate a pulse of the transverse 
waves. The lead front of the pnlse represents a whistler-like wave in upstream and 
downstream directions. 

At later time the lead front of the electromagnetic pnlse intersects the boundary of 
the computational domain, and the low-freqnency Alfven waves, which are directed by 
the external magnetic field, were observed near the spacecraft. The excitation of these 
waves are accompanied by the formation of the double current closure attached to the 
heat shield. The maximum value of the perturbations is located between the double 
currents. 

The amplitude of these waves are about of 0.05 E 0 at the linear stage and (1 — 5) E 0 


at the saturation stage. 
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These waves were generated by the gradients in plasma parameters and the Hall 
term plays the key role to support these waves. Note that such type of waves was also 
observed near a tether system moving across the magnetic field in the ionosphere, in the 
plasma environment of small moons (e.g. Io - Alfven wing), near plasma clouds moving 
across the magnetic field. 

Near the bus, the jumps in the normal component of the electric field (polarization 
field) have approximately the same value in all cases: (a) E x = 250 E 0 , E y = 50 E 0l 
E z = 60 E 0 ; (b) E x = 200 E 0l E y = 60i?o> E z = 40 E 0 ; (c) E x = 250 E 0 , E y = 40 E 0 , 
E z = 40 Eq, respectively. The tangential components of the electric field near the 
spacecraft evolve in time and their values may be different in these regimes. Note that 
the results of simulation at the saturation stage correspond to the different moments 
of time for cases (a), (b) and (c). The observed strong wave amplitudes upstream and 
downstream the SPP do not occur in the 3D hybrid simulation with a lower value of 
the /3 e = 0.000125 and /3 e = 0.00125 (see Lipatov, Sittler and Hartle, 2009). That fact 
allows us to make a conclusion that the growth of these waves is controlled by electron 
temperature and we need a further study of the mechanism of an excitation of these 
waves. 

3.2. Parallel Interaction of the Solar Wind with Solar Probe Plus. 2.5D 
simulations 

3.2.1. Parallel Interaction of the Solar Wind with Solar Probe Plus. 


2.5D simulation with large grid spacing. Since a 3D hybrid simulation of the 
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Solar Probe Plus spacecraft’s environment is very costly we have also performed a 
set of 2.5D simulations. The simulation of the oblique and quasi- parallel interaction 
shows the unphysical growth of the magnetic field barrier since the magnetic field line 
cannot penetrate though the cylindrical obstacle, ffere we discuss the results of 2.5D 
simulations with the particle absorption boundary condition (case d) and the particle 
reflection boundary condition on the heat shield (case e). 

We are also checking the results in case (d) and (e) with simulation using smaller 
grid spacing (case d’ and case e’). ffere we present only the cases with a parallel 
interaction, 9 = 0°. We use 880 x 10 6 macro-particles and a mesh with 401 x 601 grid 
points so that we have 3675 particles per cell in our 2.5D simulation. 

The simulation with reduced 2.5D configuration shows formation of an extended 
plasma wake at much longer time interval. Let us consider now the simulation results 
with j3 e = 0.1 at time, t = 2.2 T ce (t = 0.26 T transit ), where T transit is the time for particle 
transit from the left boundary to the right boundary of the computational domain. 
Figure 55 shows a proton density in the y — x plane in cases d (top, particle absorption 
boundary condition) and e (bottom, particle reflection boundary condition on the heat 
shield). We see the formation of the large void behind the heat shield and the bus. In 
case e, we also see a formation of a higher plasma density region at the front of the heat 
shield. The value of the density in this case, Fig. 55 (bottom) may exceed by factor of 2 
the value of the density in upstream in case d, Fig. 55 (top). The perturbation in the 
electromagnetic field reaches the saturation level, since the density distribution has not 


changed during the simulation. 



Figures 56, 57 and 58 show the distributions for the electric field at t = 2.2 T ce . As 
before, the perturbations in the electromagnetic field reached the saturation levels. 

The electric field upstream is calculated by the use of the simplified generalized 
Ohm ’law 

E = — (J; x B) + — (J x B), (27) 

en e c en e c 

where we dropped the VP e term because there are no strong perturbations in the 
plasma density upstream. Since the simulation does not show any perturbations in bulk 
velocity upstream we can conclude that the perturbations in the electric field E y and E z 
have a whistler-like type and these perturbations with E y ^ E z are due to the second 
term in the right side of Eq. 27 (Hall term). 

Figures 59, 60 and 61 shows 1-D cuts for the electric field profile along the x and y 
axes through the point y — 0 and x = 1.5, respectively. 

In case d, far from the spacecraft, the values of perturbation in the electric field 
components, E x , E y and E z are about of 5E X ps 5E y « ±2.0 Eq and 5E Z ps 0.1 Eq in 
the x- direction (upstream and downstream). The values of perturbation in the electric 
held components, E Xl E y and E z , are about of 5E X & 8E y ~ ±(1.0 — 2.0) E 0 (oblique 
whistler/Alfven waves or ’’shot noise”) in the y- direction, and 8E Z ps 0.02 Eq in the y- 
direction. 

In case d, near spacecraft, the E x prohle has a jump, 8E X « (±14.0, —33.0) E 0 in 
thex - direction behind the heat shield and the bus, respectively. We also see a jump in 
E y component of the electric held, SE y ~ (±15.0, —38.0) E 0 in x - direction behind the 



29 


bus and 5E y « ±55.0 E$ in y - direction, also near bottom of bus. The component E z 
has jumps, 8E Z ~ (—0.4, 0.3) E 0 in the x - direction, and 5E Z ~ —(0.4, 0.05) E 0 in the y 
- direction, Figs. 60 and 61, respectively. 

In case e, far from the spacecraft, the values of perturbation in the electric field 
components, E x , E y and E z are about of SE X ps SE y ~ ±2.0 E 0 and SE Z ~ (0.1 — 0.3) E 0 
(oblique whistlcr/Alfven waves or ’’shot noise”) in the x- direction (upstream and 
downstream), Fig. 60. The values of perturbation in the electric field components, E x , 
E y and E z , are about of 8E X « 8E y « ±1.0 E 0 (oblique whistler/Alfven waves or ’’shot 
noise”) in the y- direction and 5E Z & 0.3 Eq in the y- direction. 

In case e, near the spacecraft, the E x profile has a jump, 5E X « 500.0 E 0 and 
SE X ~ — 20.0i?o in x - direction at the front of the heat shield and behind the bus. 
We also see a jump in the E y component of the electric field, 5E y « 10.0 Eq in the x - 
direction behind the bus and 8E y « 60.0 E 0 in the y - direction. The component E z 
has jumps, 5E Z « — 1.0 Fio in the x - direction, and 8E Z « ±0.5 Eq in the y - direction, 
Figs. 60 and 61. 

The above figures demonstrate a strong oscillation in the electric field. If we assume 
the ’’shot” noise Sri/n = 0.016 then a fluctuation in the polarization electric field may 
be estimated as 5E po i = —\7p e /en e ~ 40 E 0 . This value represents the upper limit for 
the polarization electric field. 

Figures 62, 63 and 64 demonstrate distributions for the magnetic field. At the 
front of the heat shield a formation of the magnetic field barrier or build up was 
observed. At the side parts of the computational domain, compression waves were 
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observed in simulations. The value of the perturbation in the magnetic held is about of 
SB « 0.02 B 0 . 

3.2.2. Parallel Interaction of the Solar Wind with Solar Probe Plus. 
2.5D simulation with smaller grid spacing. In order to check the results of 
simulation we have repeated the simulation with a 2x smaller computational domain 
and a 2x smaller grid cell size. Let us consider now the result of simulation with fj e — 0.1 
at time, t = 1.0 T ce (t = 0.12 T transit ), where T transit is the time for particle transit 
from the left boundary to the right boundary of the computational domain. Figure 65 
shows a proton density in the y — x plane in cases d’ (top, particle absorption boundary 
condition) and e’ (bottom, particle reflection boundary condition on the heat shield). 

We see the formation of the large void behind the heat shield and the bus. As in 
case e, we see the formation of high density region at the front of the heat shield, Fig. 65 
(case e’, bottom). 

The perturbation in the electromagnetic held reaches the saturation level, since the 
density distribution has not changed during the simulation. 

Figures 66, 67, and 68 show the distributions for the electric held at t = 1.0 T ce . 
Since the density distribution had not changed during the simulation the perturbation 
in the electromagnetic held reached the saturation levels. Figures 69, 70 and 71 shows 
1-D cuts for the electric held profile along the x and y axes through the point y = 0 and 
x = 1.5 L, respectively. 

In case d’, far from the spacecraft, the values of perturbation in the electric 
held components, E x , E y and E z , are about of SE X ~ SE y ~ ±(2.0 — 3.0) E 0 
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(■ whistler/ Alfven waves or ’’shot noise”) in the x- direction (upstream and downstream), 
and SE Z ps ±(0.1 — 0.5) Eq in the x- direction. The values of perturbation in the electric 
held components, E x , E y and E~, are about of SE X & 5E y & ±(2.0 — 3.0) E 0 and 
SE Z ps ±0.05 Eq in the y- direction. 

In case d’, near spacecraft, the E x prohle has a jump, SE X ps 80.0 E 0 and 
SE X ps — 15.0i?o in x - direction at the front of the heat shield and behind the bus, 
respectively. We also see a jump in E y component of the electric held, SE y ~ 2.0 E 0 in 
x - direction behind the bus and SE y ps 100.0 E 0 in y - direction, also behind the bus. 
The component E z has jumps, SE Z ps ±2.5 Eq in x - direction in the plasma wake, and 
SE Z ps ±0.25 Eq in y - direction, Figs. 70, 71. 

In case e’, far from the spacecraft, the values of perturbation in the electric 
held components, E x , E y and E z are about of 5E X ps SE y ps ±(2.0 — 3.0) E 0 
(whistler/ Alfven waves or ’’shot noise”) in the x- direction (upstream and downstream), 
and SE Z ps ±0.01 Eq in the x- direction. The values of perturbation in the electric 
held components, E x , E y and E~, are about of 5E X & SE y ps ±(2.0 — 3.0) E 0 and 
5E Z ps ±0.1 E 0 in the y- direction. 

In case e’, near spacecraft, the E x prohle has a jump, SE X ps 150.0 Eq and 
5E X ps — 20.0i7o in x - direction at the front of the heat shield and behind the bus, 
respectively. While, 5E X ps ±20.0 E 0 in y - direction. We also see a jump in E y 
component of the electric held, 5E y ps 1.0 E 0 in x - direction behind the bus and 
SE y ps ±100.0 E 0 in y - direction also behind the bus. The component E z has jumps, 
SE Z ps 0.04 Eq in x - direction, and SE Z ps (—5.0, 1.5) E 0 in y - direction, Figs. 70, 71. 
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The above figures demonstrate a strong oscillation in the electric field. If we assume 
the ’’shot” noise Sn/n = 0.016 then a fluctuation in the polarization electric field may 
be estimated as 8E po i = —Vp e j en e « 40 E 0 . This value represents the upper limit for 
the polarization electric field. 

Figures 72, 73 and 74 demonstrate the distributions for the magnetic field. At 
the front of the heat shield a formation of the magnetic field barrier or build up was 
observed. At the side parts of the computational domain, compression waves were 
observed in simulations. The value of perturbation in the magnetic field is about of 
5B « 0.02 Bq. 

Unfortunately, 2.5D simulation was limited by the parallel interaction case. The 
simulation of the quasi-parallel and oblique interactions show the formation of the 
non-physical magnetic barrier since the magnetic filed lines cannot move through the 
cylindrical obstacle. 

Summary 2. 3D and 2.5D hybrid simulations well enough describe the global 
plasma dynamics near the SP+ spacecraft. In 3D simulation the heat shield and 
bus generate as pnlse of the transverse waves. The lead front of the pulse represents 
a whistler-like waves in upstream and downstream. In 2.5D simulations shows the 
compression-like perturbations in E field, possible associated with ’’shot” noise. The 
parallel 2.5D simulation does not demonstrate any two-current system upstream of the 
heat shield in case of the parallel interaction when the angle between the bulk velocity 
and the background magnetic field, 0u o ,b o equals zero. 

Table 3 accumulates the results of the simulations presented in this Report. Let 
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us compare the level of fluctuations in the E-helcl in 3D and 2.5D simulations. As was 
indicated on page 24 and page 28 the maximum level of fluctuations in the E-field dne 
to ’’shot” noise effects are 200 Uq x B$ in 3D simulations (case a), 100 Do x B$ (cases b 
and c), 40 Do x Bq in 2.5D simulations for the SW parameters at 9.5 R s . 

The relatively short 3D simulations produce the following oscillations in the E-field: 
5E « 5£ 0 (case a); (0. 5-1.0) E 0 (case b); 2.0 E 0 (case c), depending on the smoothing 
procedures in the E po i held or N p . In the 2.5D much longer simulations the limit for 
fluctuations in the E-field dne to a ’’shot” noise is about of 40 D 0 x B 0 . Fortunately, the 
fluctuations in the E-field in our relatively short 2.5D simulations are about of (2 — 4) Eq 
(Figs. 60, 61, 70 and 71. 

So, the single way to reduce these fluctuations is to increase the number of 
macro-particles. We also need to study the convergence of onr hybrid algorithm to 
the analytical solution in case of a strong ’’shot” noise. Possibly, we have to repeat 
onr computations with the use the advection equation applied directly to magnetic 
induction held. The scheme’s viscosity may to suppress spurious oscillations, probably 
caused by a ’’shot” noise” (see e.g. Lipatov, 2002). Note, that spurious oscillations may 
occur because our numerical model includes the second order approximation in space, 
whereas the convective term in the induction equation was presented by the hrst order 
derivative. This inconsistency may result in a spurious oscillation occurrence. 

The jump in the polarization held (2.5D simulations), E po i near the surface of SPP 
is very strong and it value does not connect with ’’shot” noise (Figs. 60, 61, 70 and 71. 

Note, that there is also possibility for an excitation of the plasma wave (c o ~ a ) pe ) 
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in the regions of the plasma wake with a low ion density, however, the study of these 
effects needs full kinetic model. 

4. Conclusions and possible future simulations 

3D hybrid simulations of the interaction of the solar wind plasma with Solar Probe 
Plus at the distance r = 9.5 Rg, have demonstrated several new features: 

• In the quasi-parallel case the current closure near a spacecraft is very complicated 
and is directed by the external magnetic field. At the front of the heat shield the 
magnetic held barrier forms, whereas strong whistler/Alfven waves occur in both 
upstream and downstream regions. At the side parts of the computational domain, 
compression waves are found. The values of the electric held oscillation near the 
spacecraft bus may be the same order as the maximum of expected electric held 
(d.c.) in an antenna (see, Table 1 from Murphy et al, 2008). 

• 3D hybrid simulations do not demonstrate the formation of the extended plasma 
void wake near the SP+ spacecraft, due to short time simulation (5 months), and 
also due to the strong polarization electric held and a large thermal ion speed 
compared with the bulk velocity. 

• For absorption boundary conditions we do not observe any perturbations in the 
plasma density upstream the SP+. 

• Now it is possible to use more than 1024 processors on ”SGI-Columbia” , 
”SGI-Plciades” (NASA Ames Center) and ’’Linux cluster Discover” (NASA 
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GSFC) parallel computer systems with hybrid MPI/OpenMP environment to 
achieve improved scalability for longer computations. 

• Simulated electric held perturbations are comparable to or exceed the maximum 
electric held expected for the SP+ spacecraft. Therefore, a wake deployment of an 
electric held plasma-wave antenna may not be a viable option for SP+. 

2.5D hybrid simulations of the interaction of the solar wind plasma with Solar 
Probe Plus at the distance r = 9.5 Rg, have demonstrated several following features: 

• 2.5D simulations show the formation of the extended plasma wake with a small 
density cavern around the spacecraft. The large spatial wake results in the 
reduction in the polarization held near the spacecraft and behind the spacecraft 
bus. 

• 2.5D simulations with a particle reflection boundary condition at the heat shield 
show the formation of the enhanced density in the external region of the plasma 
wake. Upstream of the SP+, the plasma velocity distribution may change a lot 
due to possible two-stream (beam-beam) instability. 

Future work that is still needed: 

Hybrid simulations 

• Simulations with low value of Mach number and various boundary conditions 
on the spacecraft bus and trusses conhguration; in, particular, simulations with 
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reflected boundary condition for particle at the surface of the spacecraft in order 
to take into account the effect of the spacecraft charging. 

• Simulation with much higher space resolution and a number of macro-particles to 
reduce a ’’shot” noise in density and polarization electric field; 

• 2.5D simulations with direct numerical time integration of the equation for 
magnetic field (numerical scheme for advection equation with an alternative 
direction implicit (ADI) method). Such type of scheme was used for 2.5 simulation 
by Lipatov for the study of the interaction plasma flow with a comet, Venus etc 
(see e.g., Lipatov 2002). It may possibly help us to reduce the short wavelength 
electromagnetic oscillation by the scheme resistivity. Currently, we use an implicit 
time integration scheme suggested in 1992 (limited number copies of preprint from 
the LLNL) (see e.g., Anderson and Shumaker, 1995 and further modifications in 
Lipatov, 2002). Another possibility is to use the full kinetic code that provides the 
absorption of the whistler-like wave by the wave-electrons interactions. 

• Simulations with using the (implicit or semi-implicit) equation for the time 
evolution of the electron pressure. This way may result to reduction of the 
oscillations in the electric field. 

• 2.5D simulations with fine spatial resolution to confirm the results produced in 3D 
simulations. 

• 3D hybrid long time simulation with fine resolution and the best boundary 
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conditions that minimize electromagnetic perturbations near the spacecraft bus. 

The results of these simulations will be included in future publications. 

Full kinetic simulations 

• Simulations with a test-particle approximation for electrons while the electro- 
magnetic fields will be taken from 3D hybrid simulations. Such computations 
will provide electron dynamics near the spacecraft, and provide estimates of the 
particle flux (velocity distribution) at the surface of the spacecraft. 

• Simulation that will take into account the charging of the spacecraft, charge 
separation effects, outgassing from the heat shield, and effects due to 
photoionization and electron impact ionization effects near the spacecraft. 

Hall-MHD or two-fluid simulations 

• Hall-MHD or two-fluid simulations. These models have no any ’’shot 1 ' noise effects 
and may produce a perfect current closure and the low-frequency wave system 
near the SPP. 

• Two-fluid models (separate fluid approximation for ions and electrons) may help 
us to study the low- and high-frequency plasma oscillations inside the plasma 
wake in the regions with a low value of the density. 
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Table 1 . Important values of the electromagnetic field in an antenna from 
” Thermal Design Consideration for the Solar Probe Electric Field Antenna 
by Murphy et ah, Feb. 2008”. 



4 R s 

12i? s 

20 R s 

U xB 

12 V/m 

2.5 V/m 

0.5 V/m 

max i DC 

lV/m 

0.1 V/m 

0.1 V/m 

Emax , 100kHz 

10“ 3 V/m Hz" 1 / 4 

10" 4 V/m Hz" 1 / 4 

10" 4 V/m Hz" 1 / 4 

V Noise at 100 kHz 

10“ 6 V/Hz- 1 / 4 

3 x 10 -7 V/Hz -1 / 4 

2 x 10 -7 V /Hz -1 / 4 
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Table 2. Solar Wind Parameters for Modeling. 


Case Uq 

nsw 

B 0 

@bu 

Ma P p 

Pe 

Uo x Bo 

N p per cell 

Comments 

a 

200 krn/s 

5 x 10 3 

1500 nT 

n° 

1.5 

0.1 

0.1 

0.3 V/m 

75 

no smoothing in E po i 

b 

200 krn/s 

5 x 10 3 

1500 nT 

n° 

1.5 

0.1 

0.1 

0.3 V/m 

800 

smoothing in E po i 

c 

200 km/s 

5 x 10 3 

1500 nT 

n° 

1.5 

0.1 

0.1 

0.3 V/m 

800 

no smoothing in E po i 

d 

200 km/s 

5 x 10 3 

1500 nT 

0° 

1.5 

0.1 

0.1 

0.3 V/m 

800 

2.5D, with absorption 

d’ 

200 km/s 

5 x 10 3 

1500 nT 

0° 

1.5 

0.1 

0.1 

0.3 V/m 

800 

2.5D, with absorption; 











higher spatial resolution 

e 

200 km/s 

5 x 10 3 

1500 nT 

0° 

1.5 

0.1 

0.1 

0.3 V/m 

800 

2.5D with reflection 

e’ 

200 km/s 

5 x 10 3 

1500 nT 

0° 

1.5 

0.1 

0.1 

0.3 V/m 

800 

2.5D with reflection; 


higher spatial resolution 












Case t(T ce ) 

SE X 

5Ey 

5E Z 

SB X 

5By 

5B Z 

E n at bus 

a 

0.014 

0 . 02 - 0.15 

0 . 03 - 0.1 

0 . 02 - 0.1 

0.001 

0.0002 

0.001 

0 . 01 - 0.1 

a 

0.29 

5.0 

5.0 

5.0 

0.005 

0.01 

0.1 

60-250 

b 

0.014 

0.2 

0.2 

0.2 

0.002 

0.005 

0.02 

0 . 4 - 1.0 

b 

0.15 

0 . 5 - 1.0 

0 . 5 - 1.0 

0 . 5 - 1.0 

0.005 

0.01 

0.1 

40-200 

c 

0.29 

2.0 

2.0 

2.0 

0.005 

0.01 

0.1 

35-250 

d 

2.2 

1 . 0 - 2.0 

1 . 0 - 2.0 

0.1 

0.02 

0.02 

0.005 

30-60 

d ’ 

2.2 

2 . 0 - 3.0 

2 . 0 - 3.0 

0.1 

0.02 

0.02 

0.005 

30-60 

e 

2.2 

1 . 0 - 2.0 

1 . 0 - 2.0 

0.3 

0.02 

0.02 

0.005 

30-60 

e ’ 

2.2 

2 . 0 - 3.0 

2 . 0 - 3.0 

0.3 

0.02 

0.02 

0.005 

30-60 















Figure 1 . Scheme of the interaction of the solar wind with SP+. The spiral magnetic 
field is inside the x-z plane. (1) a heat shield; (2) a support system; (3) a spacecraft bus. 
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Figure 2. Solar wind ion density in the y-x (z = 0) and z-x (y = 0) planes. U 0 = 
200km/s, Ma = 1.5, B 0 = 1500 nT, E 0 = UqBq = 0.3 V/m, nsw = 5 x 10 3 cm~ 3 , 
/3 P = 0.1, (3 e = 0.1, 9bu = 11°. Linear perturbations occur at the beginning of simulation, 


at t — 0.014 T ce , (case a). 
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Figure 3. Solar wind ion density in the y-x (z = 0) and z-x (y = 0) planes. U 0 = 
200 km/s, M a = 1.5, B 0 = 1500 nT, n = 5 x 10 3 cm- 3 , = 0.1, (3 e = 0.1, 6 bu = 11°. 

Linear perturbations occur at the beginning of simulation, at t = 0.014 T ce , (case b). 
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Figure 4. Electric field component E X (E 0 ) in the y-x ( z = 0) and z-x (y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, E 0 = UqBq = 0.3V/m, nsw = 5 x 1() 3 crrf 3 , 
/3 P = 0.1, f3 e = 0.1, dbu = 11°. Linear perturbations at the beginning of simulation, at 


t = 0.014 T ce , (case a). 
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Figure 5. Electric field component E X (E 0 ) in the y-x (z = 0) and z-x (r/ = 0) planes. 
Eo = 200km/s, Myi = 1.5, B 0 = 1500nT, E 0 = UqBq = 0.3V/m, nsw = 5 x 1() 3 crrf 3 , 
/3 P = 0.1, /3 e = 0.1, dbu = 11°. Linear perturbations at the beginning of simulation, at 
t = 0.014 T ce (case b). 
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Figure 6. Electric field component E y (E 0 ) in the y-x (z = 0) and z-x (y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm~ 3 , (3 P = 0.1, j3 e = 0.1, 9j, u = 11°. 
Linear perturbations at the beginning of simulation, at t — 0.014 T ce , (case a). 
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Figure 7. Electric field component E y (E 0 ) in the y-x (z = 0) and z-x (y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm -3 , (3 P = 0.1, (3 e = 0.1, 9j, u = 11°. 
Linear perturbations at the beginning of simulation, at t — 0.014 T ce (case b). 
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Figure 8. Electric field component E Z (E 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm~ 3 , (3 P = 0.1, (3 e = 0.1, 9j, u = 11°. 
Linear perturbations at the beginning of simulation, at t — 0.014 T ce , (case a). 
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Figure 9. Electric field component E Z (E 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm~ 3 , (3 P = 0.1, (3 e = 0.1, 9j, u = 11°. 
Linear perturbations at the beginning of simulation, at t — 0.014 T ce (case b). 
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Figure 10. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (y = 0, 
z = 0). Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , /3 p = 0.1, (3 e = 0.1, 
9b u = 11°. Linear perturbations at the beginning of simulation, at t — 0.014 T ce , (case a). 
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Figure 11. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) ( y = 0, 
z = 0). U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , /3 p = 0.1, (3 e = 0.1, 
9b u = 11°. Linear perturbations at the beginning of simulation, at t = 0.014 T ce (case b). 
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Figure 13. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (x = 
1.5, z = 0). U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , (3 V = 0.1, 
(3 e = 0.1, 9bu = 11°. Linear perturbations at the beginning of simulation, at t = 0.014 T ce 
(case b). 
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Figure 14. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (x = 
1.5, y = 0). U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , fi p = 0.1, 


(3 e = 0.1, 0f- }U — 11°. Linear perturbations at the beginning of simulation, at t — 0.014 T ce , 
(case a). 
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Figure 15. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (x = 
1.5, y = 0). U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , (3 V = 0.1, 
/3 e = 0.1, 9b u = 11°. Linear perturbations at the beginning of simulation, at t = 0.014 T ce 
(case b). 
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Figure 16. Magnetic field component B X (B 0 ) in the y-x (z = 0) and z-x (y = 0) planes. 
U 0 = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm” 3 , /3 P = 0.1, /3 e = 0.1, 9b u = 11°. 
Linear perturbations at the beginning of simulation, at t = 0.014 T ce , (case a). 
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Figure 17. Magnetic field component B X (B 0 ) in the y-x (z = 0) and z-x (y = 0) planes. 
U 0 = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm” 3 , /3 P = 0.1, /3 e = 0.1, 9b u = 11°. 
Linear perturbations at the beginning of simulation, at t = 0.014 T ce (case b). 
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Figure 18. Magnetic field component B y (B 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n = 5xl0 3 cm -3 , /3 P = 0.1, /3 e = 0.1, 9b u = 11°. 
Linear perturbations at the beginning of simulation, at t = 0.014 T ce , (case a). 
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Figure 19. Magnetic field component B y (B 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n = 5xl0 3 cm -3 , /3 P = 0.1, /3 e = 0.1, 9b u = 11°. 
Linear perturbations at the beginning of simulation, at t = 0.014 T ce (case b). 
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Figure 20. Magnetic field component B Z (B 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm -3 , /3 P = 0.1, /3 e = 0.1, 6f m = 11°. 
Linear perturbations at the beginning of simulation, at t = 0.014 T ce , (case a). 
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Figure 21. Magnetic field component B Z (B 0 ) in the y-x (z = 0) and z-x (?/ = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm” 3 , /3 P = 0.1, /3 e = 0.1, 6f m = 11°. 
Linear perturbations at the beginning of simulation, at t = 0.014 T ce (case b). 
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Figure 22. Solar wind ion density in the y-x (z = 0) and z-x ( y 


0) planes. 


U 0 = 200 krn/s, Ma = 1.5, B$ = 1500nT, E 0 = UqB 0 = 0.3V/m, risw = 5 x 10 3 cm 3 , 
/ 3 P = 0.1, /3 e = 0.1, db u = 11°. Nonlinear saturation of the perturbations at t — 


0 . 1 04 T ce (0.013 T transit ) , (case a). 
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Figure 23. Solar wind ion density in the y-x (z = 0) and z-x ( y = 0) planes. 
U 0 = 200km/s, Ma = 1.5, B 0 = 1500nT, E 0 = U 0 B 0 = 0.3V/m, nsw = 5 x 10 3 cm -3 , 
/3 P = 0.1, /3 e = 0.1, 9b u = 11°. Nonlinear saturation of the perturbations at t = 


0.15 T ce (0.02 T transit ) (case b). 
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Figure 24. Solar wind ion density in the y-x (z = 0) and z-x ( y = 0) planes. 
U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, E 0 = U 0 B 0 = 0.3V/m, nsw = 5 x 10 3 cm -3 , 
/3 P = 0.1, /3 e = 0.1, 9bu = 11°. Nonlinear saturation of the perturbations at t = 


0.29 T ce (0.04 T transit ) (case c). 
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Figure 25. Electric field component E X (E 0 ) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, E 0 = UqB 0 = 0.3 V/m, nsw = 
5 x 10 3 cm~ 3 , /3 P = 0.1, /3 e = 0.1, 6b u = 11°. Nonlinear saturation of the perturbations at 


t = 0.104 T ce (0.013 T transit ) (case a). 
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Figure 26. Electric field component E X (E 0 ) in the y-x (z = 0), z-x ( y = 0) and y-z 
(x = 0) planes. U 0 = 200km/s, M 4 = 1.5, B 0 = 1500 11 T, E 0 = U 0 B 0 = 0.3 V/m, 
risw = 5 x 10 3 cm -3 , (3 P = 0.1, f3 e = 0.1, 6i >u = 11°. Nonlinear saturation of the 
perturbations at t = 0.15T ce (0.02T transit ) (case b). 
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Figure 27. Electric field component E X (E 0 ) in the y-x (z = 0) and z-x (y = 0) planes. 
U 0 = 200km/s, Ma = 1.5, nsw = 5 x 10 3 cm -3 , (3 P = 0.1, (3 e = 0.1, 9 bu = 11°. Nonlinear 
saturation of the perturbations at t — 0.29 T ce (0.04 Transit) (case c). 
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Figure 28. Electric field component E y (E 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm -3 , (3 P = 0.1, (3 e = 0.1, 0f m = 11°. 
Nonlinear saturation of the perturbations at t = 0.104 T ce (0. 013 T tra nsit) (case a). 
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Figure 29. Electric field component E y (E 0 ) in the (z = 0), z-x ( y = 0) and y-z 
(x = 0) planes. f/ 0 = 200km/s, Ma = 1.5, E> 0 = 1500 nT, u = 5x 10 3 cm~ 3 , /3 P = 0.1, 

P e = 0.1, 6b u = 11°. Nonlinear saturation of the perturbations at t — 0.15 T ce (0.02 T transit ) 

(case b). 
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Figure 30. Electric field component E y (E 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
U 0 = 200km/s, M a = 1.5, B 0 = 1500 nT, n = 5x 10 3 cm- 3 , f3 p = 0.1, (3 e = 0.1, 6 bu = 11°. 
Nonlinear saturation of the perturbations at t = 0.29 T ce (0.04 T tran sit) (case c). 
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Figure 31. Electric field component E Z (E 0 ) in the y-x (z = 0) and z-a; (r/ = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm~ 3 , (3 P = 0.1, f3 e = 0.1, 9b u = 11°. 


Nonlinear saturation of the perturbations at t = 0.104 T ce (0. 013 Transit) (case a). 
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Figure 32. Electric field component E Z (E 0 ) in the y-x (z = 0), z-x (y = 0) and y-z 
( x = 0) planes. U 0 = 200km/s, Ma = 1.5, n — 5 x 10 3 cm -3 , /3 P = 0.1, (3 e = 0.1, 
dbu = 11°. Nonlinear saturation of the perturbations at t — 0.15T ce (0.02T trans i t ) (case 
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Figure 33. Electric field component E Z (E 0 ) in the y-x (z = 0) and z-x (■ y = 0) planes. 
U 0 = 200km/s, M a = 1.5, B 0 = 1500 nT, n = 5x 10 3 cm~ 3 , f3 p = 0.1, (3 e = 0.1, 6 bu = 11°. 
Nonlinear saturation of the perturbations at t — 0.29 T ce ( 0.04 T trans it) (case c). 
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Figure 34. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (y = 0, 
z = 0). Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , /3 p = 0.1, (3 e = 0.1, 
9bu = 11°. Nonlinear saturation of the perturbations at t = 0.104 T ce (0. 013 Transit) (case 



z = 0). U 0 = 200km/s, Ma = 1.5, n — 5 x 10 3 cm 3 , (3 P = 0.1, (3 e = 0.1, 9 bu = 11°. 


Nonlinear saturation of the perturbations at t — 0.15 T ce (0.02 Transit) (case b). 
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Figure 37. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (y = 0, 
z = 0). Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , /3 p = 0.1, (3 e = 0.1, 
9bu = 11°. Nonlinear saturation of the perturbations at t = 0.104 T ce (0. 013 Transit) (case 
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Figure 38. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) ( y = 0, 
z = 0). Uq = 200km/s, Ma = 1.5, n — 5 x 10 3 cm -3 , /3 P = 0.1, (3 e = 0.1, 6 bu = 11°. 
Nonlinear saturation of the perturbations at t — 0.15 T ce (0.02 Transit) (case b). 
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Figure 39. 1-D cuts for the electric field component E x { 
z = 0). U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 
Obu = 11°. Nonlinear saturation of the perturbations at 
c). 


E 0 ), E y (E 0 ) and E Z (E 0 ) (y = 0, 
x 10 3 cm -3 , P p = 0.1, P e = 0.1, 
t = 0.29 X^e (0.04 T trans it) (case 





Figure 40. 1-D cuts for the electric field component E X (E 0 , E y (E 0 ) and E Z (E 0 ) (x 
z = 0). Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , /3 p = 0.1, (3 e 
6bu = 11°. Nonlinear saturation of the perturbations at t = 0.104 T ce (0. 013 Transit 




Figure 41. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (x 
1.5, z = 0). Uq = 200km/s, Ma = 1.5, n — 5 x 10 3 cm -3 , f3 p = 0.1, /3 e = 0.1, 9 bu = F 
Nonlinear saturation of the perturbations at t — 0.15 T ce (0.02 Transit) (case b). 





Figure 42. 1-D cuts for the electric field component E X {E 0 ), E y (E 0 ) and E Z {E 0 ) (x = 
1.5, z = 0). Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , /3 p = 0.1 
(3 e = 0.1, 6b u = 11°. Nonlinear saturation of the perturbations at t = 0. 29 T ce (0. 04 T trans it 






Figure 43. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (x = 
1.5, y = 0). Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , 0 p = 0.1, (3 e = 
0.1, 9bu = 11°. Nonlinear saturation of the perturbations at t = 0.104 T ce (0 . 013 T tran sit 





Figure 44. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (x 
1.5, y = 0). Uq = 200km/s, Ma = 1.5, n — 5 x 10 3 cm -3 , f3 p = 0.1, /3 e = 0.1, 9 bu = F 
Nonlinear saturation of the perturbations at t — 0.15 T ce (0.02 Transit) (case b). 
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Figure 45. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 and E Z (E 0 ) (x 
y = 0). Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , f3 v = 0.1, {3 e 
6bu = 11°. Nonlinear saturation of the perturbations at t = 0.29 T ce ( 0.04 T tran sit 
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Figure 46. Magnetic field component B X (B 0 ) in the y-x (z = 0) and z-x (y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm” 3 , (3 P = 0.1, (3 e = 0.1, 6f m = 11°. 
Nonlinear saturation of the perturbations at t = 0.104 T ce (0. 013 T tra nsit) (case a). 
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Figure 47. Magnetic field component B X (B 0 ) in the y-x (z = 0), z-x (y = 0) and y-z 
(x = 0) planes. U 0 = 200km/s, Ma = 1.5, n — 5 x 10 3 cm -3 , /3 p = 0.1, /3 e = 0.1, 
dbu = 11°. Nonlinear saturation of the perturbations at t — 0.15 ^(0.02 T transit ) (case 


b). 
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Figure 48. Magnetic field component B X (B 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
U 0 = 200km/s, M a = 1.5, B 0 = 1500 nT, n = 5x 10 3 cm~ 3 , f3 p = 0.1, (3 e = 0.1, 6 bu = 11°. 
Nonlinear saturation of the perturbations at t — 0.29T ce (0.0AT transit ) (case c). 
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Figure 49. Magnetic field component B y (B 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n = 5xl0 3 cm” 3 , (3 P = 0.1, (3 e = 0.1, 6f m = 11°. 
Nonlinear saturation of the perturbations at t = 0.104 T ce (0. 013 T tra nsit) (case a). 
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Figure 50. Magnetic field component B y (B 0 ) in the y-x (z = 0), z-x (y = 0) and y-z 
(x = 0) planes. U 0 = 200km/s, Ma = 1.5, n — 5 x 10 3 cm -3 , /3 P = 0.1, (3 e = 0.1, 
dbu = 11°. Nonlinear saturation of the perturbations at t — 0.15 ^(0.02 T transit ) (case 


b). 
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Figure 51. Magnetic field component B y (B 0 ) in the y-x (z = 0) and z-x ( y = 0) planes. 
U 0 = 200km/s, M a = 1.5, B 0 = 1500 nT, n = 5x 10 3 cm~ 3 , f3 p = 0.1, (3 e = 0.1, 6 bu = 11°. 
Nonlinear saturation of the perturbations at t — 0.29 T ce ( 0.04 Transit) (case c). 
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Figure 52. Magnetic field component B Z (B 0 ) in the y-x (z = 0) and z-x (?/ = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 1500nT, n = 5xl0 3 cm -3 , (3 P = 0.1, (3 e = 0.1, 6f m = 11°. 
Nonlinear saturation of the perturbations at t = 0.104 T ce (0. 013 T tra nsit) (case a). 
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Figure 53. Magnetic field component B Z (B 0 ) in the y-x (z — 0), z-x ( y = 0) and y-z 
(x = 0) planes. U 0 = 200km/s, M 4 = 1.5, n = 5 x 10 3 cm -3 , /3 P = 0.1, (3 e = 0.1, 
dbu = 11°. Nonlinear saturation of the perturbations at t = 0.15 T ce ( 0.02 Transit) (case 


b). 
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Figure 54. Magnetic field component B Z (B 0 ) in the y-x (2 = 0) and 2 -x (y = 0) planes. 
U 0 = 200km/s, M a = 1.5, B 0 = 1500 11 T, n = 5x 10 3 cm~ 3 , f3 p = 0.1, (3 e = 0.1, 6 bu = 11°. 
Nonlinear saturation of the perturbations at t = 0.29 T ce (0.04 T tran sit) (case c). 
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Figure 55. Solar wind ion density in the y-x plane. Uq = 200km/s, Ma = 1.5, 
Bq = 1500 nT, Eq = UqB 0 = 0.3 V/m, nsw = 5 x 10 3 cm -3 , f3 p = 0.1, /3 e = 0.1, 9 bu = 0°. 
Nonlinear saturation of the perturbations at t — 2.2 T ce (0.26 Transit) (case d, top; case e, 
bottom). 
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Figure 56. Electric field component E X (E 0 ) in the y-x plane. U$ = 200km/s, = 1.5, 

£>0 = 1500 nT, F 0 = U 0 B 0 = 0.3 V/m, nsw = 5 x 10 3 cm -3 , (3 P = 0.1, /3 e = 0.1, 6 bu = 0°. 
Nonlinear saturation of the perturbations at t = 2.2 T ce ( 0.26 Transit) (case d, top; case e, 
bottom). 
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Figure 57. Electric field component E y (E 0 ) in the y-x plane. Uq = 200km/s, = 

1.5, rtsw = 5 x 10 3 cm~ 3 , (5 p = 0.1, (3 e = 0.1, 9 bu = 0°. Nonlinear saturation of the 


perturbations at t — 2.2 T ce (0.26 T tra nsit) (case d, top; case e, bottom). 





Y/L Y/L 


Ez/Eo 



-10 -5 0 5 10 

X/L 


Ez/Eo 



0.10 

0.05 

0.00 

- 0.05 

- 0.10 


-10 -5 0 5 10 

X/L 


Figure 58. Electric field component E Z (E 0 ) in the y-x plane. Uq = 200km/s, M 4 = 
1.5, rtsw = 5 x 10 3 cm~ 3 , (3 p = 0.1, /3 e = 0.1, 9b u = 0°. Nonlinear saturation of the 
perturbations at t = 2.2T ce (0.26T tra nsit) (case d, top; case e, bottom). 
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Figure 59. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (y — 0). 
U 0 = 200km/s, M A = 1.5, B 0 = 1500 nT, n = 5 x 10 3 cm- 3 , p p = 0.1, p e = 0.1, 6 bu = 0°. 
Nonlinear saturation of the perturbations at t = 2.2 T ce (0.26 T transi t) (case d, left; case e, 
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Figure 60. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (y — 0). 
U 0 = 200km/s, M a = 1.5, B 0 = 1500 nT, n = 5 x 10 3 cm' 3 , (3 p = 0.1, (3 e = 0.1, 9 bu = 0°. 


Nonlinear saturation of the perturbations at t = 2.2 T ce (0.26T transit ) (case d, left; case e, 
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Figure 61. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (x = 
1.5). U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , (3 V = 0.1, f3 e = 0.1, 
Obu = 0°. Nonlinear saturation of the perturbations at t = 2.2 T ce (0.26 T tran sit) (case d, 


left; case e, right). 
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Figure 62. Magnetic field components B X (B 0 ) in the y-x plane. t/ 0 = 200km/s, = 

1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , /3 p = 0.1, /3 e = 0.1, 9^ = 0°. Nonlinear saturation 
of the perturbations at t = 2.2 T ce (0. 26 T trans it) (case d, top; case e, bottom). 
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Figure 63. Magnetic field components B y (B 0 ) in the y-x plane. U 0 = 200km/s, Ma = 
1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , (3 p = 0.1, /3 e = 0.1, 6^ = 0°. Nonlinear saturation 
of the perturbations at t = 2.2 T ce (0. 26 T trans it) (case d, top; case e, bottom). 
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Figure 64. Magnetic field component B Z (B 0 ) in the y-x plane. Uq = 200km/s, M a — 


1.5, B 0 = 1500 nT, n — 5 x 10 3 cm 3 , /3 P = 0.1, /3 e = 0.1, 9 bu = 0°. Nonlinear saturation 


of the perturbations at t = 2.2 T ce (0. 26 T trans it) (case d, top; case e, bottom). 
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Figure 65. Solar wind ion density in the y-x plane. Uq = 200km/s, Ma = 1.5, 
Bq = 1500 nT, E 0 = U 0 B 0 = 0.3 V/m, nsw = 5 x 10 3 cm -3 , (3 P = 0.1, /3 e = 0.1, 6 bu = 0°. 
Nonlinear saturation of the perturbations at t — 1.0 T ce ( 0.12 T tran sit) (case d’, top; case 
e’, bottom). 
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Figure 66. Electric field component E X (E 0 ) in the y-x plane. Uq = 200km/s, = 1.5, 

i?o = 1500 nT, F 0 = U 0 B 0 = 0.3 V/m, nsw = 5 x 10 3 cm -3 , (3 P = 0.1, /3 e = 0.1, 6 bu = 0°. 
Nonlinear saturation of the perturbations at t = 1.0 T ce ( 0.12 T tran sit) (case d’, top; case 
e’, bottom). 
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Figure 67. Electric field components E y (E 0 ) in the y-x plane. Uq = 200km/s, Ma = 
1.5, rtsw = 5 x 10 3 cm -3 , = 0.1, (3 e = 0.1, 9 bu = 0°. Nonlinear saturation of the 

perturbations at t = 1.0 T ce (0. 12 Transit) (case d’, top; case e’, bottom). 
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Figure 68. Electric field components E Z (E 0 ) in the y-x plane. Uq = 200km/s, Ma = 
1.5, rtsw = 5 x 10 3 cm~ 3 , /3 p = 0.1, (3 e = 0.1, 9 bu = 0°. Nonlinear saturation of the 
perturbations at t = 1.0 T ce (0. 12 Transit) (case d’, top; case e’, bottom). 
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Figure 69. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (y — 0). 
U 0 = 200km/s, M A = 1.5, B 0 = 1500 nT, n = 5 x 10 3 cm- 3 , p p = 0.1, p e = 0.1, 6 bu = 0°. 
Nonlinear saturation of the perturbations at t = 1.0 T ce ( 0.12 T transi t) (case d’, left; case 
ght) 
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Figure 70. 1-D cuts for the electric field con 
U 0 = 200km/s, Ma = 1.5, B 0 = 1500 nT, n 
Nonlinear saturation of the perturbations a 
ight) 
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Figure 71. 1-D cuts for the electric field component E X (E 0 ), E y (E 0 ) and E Z (E 0 ) (x 
1.5). Uq = 200km/s, Ma = 1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , /3 P = 0.1, j3 e = 0 
Obu = 0°. Nonlinear saturation of the perturbations at t — 1.0 T ce (0.12 T transit ) (case < 


left; case e’, right). 
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Figure 72. Magnetic field components B X (B 0 ) in the y-x plane. t/ 0 = 200km/s, = 

1.5, i?o = 1500 nT, n = 5 x 10 3 cm -3 , /3 p = 0.1, /3 e = 0.1, = 0°. Nonlinear saturation 


of the perturbations at t — 1 . 0 T ce (0.12 T transit ) (case d’, top; case e’, bottom). 
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Figure 73. Magnetic field components B y (B 0 ) in the y-x plane. U Q = 200km/s, Ma = 
1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , (3 p = 0.1, /3 e = 0.1, 9f, u = 0°. Nonlinear saturation 
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of the perturbations at t = 1.0 T ce ( 0.12 T trans it) (case d’, top; case e’, bottom). 
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Figure 74. Magnetic field component B Z (B 0 ) in the y-x plane. Uq = 200km/s, Ma = 
1.5, B 0 = 1500 nT, n — 5 x 10 3 cm -3 , (3 p = 0.1, /3 e = 0.1, 6^ = 0°. Nonlinear saturation 
of the perturbations at t = 1.0 T ce ( 0.12 T trans it) (case d’, top; case e’, bottom). 
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